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We obtain the complete phase diagram of the antiferromagnetic J1-J2 model, 0 < a = J2/JI < 1 , 
within the framework of the 0(N) nonlinear sigma model. We find two magnetically ordered 
phases, one with Neel order, for a < 0.4, and another with collinear order, for a > 0.6, separated 
by a nonmagnetic region, for 0.4 < a < 0.6, where a gapped spin liquid is found. The transition at 
a = 0.4 is of the second order while the one at a = 0.6 is of the first order and the spin gaps cross at 
a = 0.5. Our results are exact at N —)> 00 and agree with numerical results from different methods. 
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Quantum phase transitions (QPTs) occur when the 
ground state properties of a certain physical system un¬ 
dergo dramatic changes as one, or more, internal or ex¬ 
ternal, parameters are varied [T]. Examples include, but 
are not restricted to, magnetic phase transitions between 
two distinct magnetic ground states or between a mag¬ 
netic state and a nonmagnetic one, driven for example 
by an applied field, pressure, or the coupling to oder de¬ 
grees of freedom. QPTs are usually labelled according to 
the behaviour of some order parameter (OP) close to the 
quantum critical point (QCP) |2], and are said to be of 
the second order (2nd order) when the OP vanishes con¬ 
tinuously as the QCP is approached, or of the first order 
(1st order) when the OP has a finite value near the QCP 
and jumps discontinuously to zero above it. Furthermore, 
knowledge of the range of the interactions, symmetries 
of the Hamiltonian and dimension of the OP, allow us to 
classify QPTs into universality classes [2], and help us to 
wirte down a Landau-Ginzburg free energy (LGFE) to 
describe such phase transitions (PTs). Typically, LGFEs 
up to the 4th power of the OP are enough to describe a 
2nd order PT, while LGFEs up to the 6th power of the 
OP are necessary to describe a 1st order PT. 

The 0(N ) quantum nonlinear sigma model (NLSM) 
has long been acknowledged to be a very convenient 
framework to describe 2nd order magnetic PTs in spin 
systems, such as, for example, the antiferromagnetic 
(AF) Heisenberg Hamiltonian, in two dimensions, with 
nearest-neighbour interactions on a square lattice [3]. 
Here the QPT occurs between a Neel ordered magnetic 
ground state, where the OP is the sublattice magnetiza¬ 
tion, a / 0, and a nonmagnetic state (cr = 0) with a finite 
spin gap, A / 0, as the OP at zero temperature. Such 
transition is driven by quantum fluctuations set by some 
coupling constant, g , and is of the 2nd order, as both a 
and A vanish continuously at the QCP, g c . Despite being 
nonlinear, at the mean field level (N -A 00 ) the model is 
quadratic, exactly solvable, and produces the usual mean 
field values for the critical exponents of the Heisenberg 
universality class, a oc [g c — g)P, for the ordered regime 
{9 < 9c ), with P = 1/2, and A oc {g — g c ) u > f° r 
nonmagnetic phase (g > g c ), with v = 1 [4]. 


First order PTs in spin systems occur whenever two 
magnetic phases cannot be continuously connected to one 
another by some order parameter. This is what happens, 
for example, already at the classical level, between the 
Neel- and collinear-type ordering phases of the J\ — J 2 
model, at the border a = 0.5. When quantum fluctua¬ 
tions are taken into account, a nonmagnetic region opens 
up around a = 0.5 [5J and a gapped spin liquid phase is 
found for 0.4 < a < 0.6 [6j. Although the precise nature 
of the nonmagnetic state is still under debate (typical 
candidates range from dimer to plaquette or VBS phases) 
the nature of the transition at a = 0.4 is agreed to be of 
the 2nd order by either numerical and theoretical meth¬ 
ods, like for example the NLSM |7j. For the transition at 
a = 0.6, different numerical techniques, including series 
expansion (8|, quantum Monte Carlo |9j. exact diagonal- 
ization [lU], and DMRG [TT] , strongly indicate it to be of 
the 1st order m, but from the theoretical point of view 
no conclusive statement has yet been presented. More 
importantly, this poses serious questions on the applica¬ 
bility of the NLSM to describe a 1st order PT in frus¬ 
trated magnetic systems m , specially since no unusual 
powers of the OP are to be expected. 

In this work we derive and apply the O(N) NLSM for¬ 
malism for the J\ — J 2 Heisenberg model, for the whole 
range of parameters, 0 < a < 1. Up to the classi¬ 
cal border, 0 < a < 0.5, the model describes smooth 
fluctuations of the staggered order parameter on top of 
a Neel ordered ground state and possesses a 2nd order 
phase transition, at a = 0.4, driven by quantum fluctua¬ 
tions, towards a nonmagnetic, gapped spin liquid phase. 
Beyond the classical border, 0.5 < a < 1, the model 
describes, instead, smooth fluctuations of the staggered 
order parameter on top of a collinearly ordered ground 
state. Remarkably, although at the mean field (N -A 00 ) 
level the model remains quadratic and exactly solvable, 
we show that its quantum dynamics is importantly mod¬ 
ified by a term proportional to the AF order parameter, 
which causes significant changes on the behaviour of the 
OP at zero temperature. The nonmagnetic, gapped spin 
liquid and collinear phases can no longer be continuously 
connected and a 1st order QPT is theoretically obtained. 
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FIG. 1: Collinear state of the J\ — J 2 model for 0.6 < a < 1 
as two interpenetrated Neel structures A and B , with nearest 
neighbour, Ji, and next-to-nearest neighbour, J 2 , AF inter¬ 
actions, upon which we shall build up our NLSM. 


The J\ — J 2 Heisenberg spin-Hamiltonian is given by 

H = -hY, S*' 7 + -h ^Sr Sj, (1) 

(iJ) « iJ » 

where Ji > 0 and J 2 > 0 are, respectively, the AF su¬ 
perexchanges between nearest-neighbors, and next- 

to-nearest neighbors, ((i, j)), spins on a two dimen¬ 
sional square lattice. The Hamiltonian 0 exhibits two 
types of magnetic order: Neel order, with wave vector at 
q = (7r,7r), for a < 0.4, and collinear order, with wave 
vectors at q = (7r, 0) and/or q = (0,7r), for a > 0.6 [14]. 

For the Neel ( 7 r, 7r) phase, different effective field theo¬ 
ries, of the NLSM type, have been proposed [7], and they 
all succeed in describing the 2nd order PT at a = 0.4. 
For the collinear ( 77 , 0) and/or (0,7r) phase, instead, no 
such description has been provided yet, and we shall pro¬ 
ceed as follows: we treat the collinear magnetic state 
as a result of two interpenetrated Neel ordered sublat¬ 
tices and introduce an double coherent spin-state ba¬ 
sis, with spin operators labeled by indices A and H, 
see Fig. [l] We then associate the spins operators in 
Eq. 0 to vector fields n A {i) and %(i) that describe 
long wavelength deviations from the Neel state in each 
sublattice. As usual we parametrize the spin-1 fields 
into a smooth, m, and a fast and uniform, L, varying 
components, n A ,B = ^a,b^a,b \/1 - {clL a ,b) 2 + uL a ,b, 
where 0 A ,B(i) = +1 for spin f and = —1 for spin 
and a = a d /S. To satisfy n A: B • n A ,B = 1 we as¬ 
sume that TTi A B • ™a,b ~ 1, while fh A ^B * L a ,b ~ 0, 
o?L a ,b • L A: b « 1, and rh A • Lb = itlb • L A = 0. After 


integration over L the action for the smooth fields is 


5 = ^ 


/ 


Ps 
2 

7o (rh B 


(X7fh A ) 2 + (Vm B ) 2 + 


(d r fh A ) 2 + ( d r fh B ) 2 


d x d y rh A + m A • d x d y m B ) 

- c{ 2 [m A x 5 r mA) * {m B x 9,%)] , ( 2 ) 

where = 2J 2 5 2 is the spin stiffness in two dimen¬ 
sions, Co = \/2Sa^/l6J2 — J 2 and c\ — are spin- 

wave velocities, and 70 = yG The first line in Eq. (p| 
corresponds to the usual NLSM for the two Neel suo- 
structures of Fig. [l] labelled A and B, which are decou¬ 
pled when J\ = 0. For J\ 7 ^ 0, however, two couplings 
arise: the first one involves only gradient terms and pro¬ 
duces different spin-wave velocities along the diagonals 
m ; the second, and more important one, is a result of 
the coupled pressession of magnetic moments on the two 
Neel sub-structures and modifies importantly the dynam¬ 
ics of the problem, ultimately leading to the first order 
character of the phase transition at a = 0 . 6 . 

In the magnetically ordered phase we can write 
friA.B = Tc,(a,b)£ + 'R‘y,(A,B)y + The 7 r fields are 
associated to the quantum fluctuations and the a field to 
the staggered OP. We introduce the Lagrange multiplier 
Svinc oc j iX(\m A \ — l)+iA(|ms| —1) and after integrat¬ 
ing out transverse fluctuations we end up with the parti¬ 
tion function Z{f3) = A/ 7 / D[iX]D[a]e~ NSe ^^ x,cr \ where 


^eff [^5 — 


N ■ 


N 


—Tr In (A(<9) + iXI) + J — iX(a 2 


geo 


is given in terms of 


m = 


ai(d) 0 Ci 2 {d) 0 

0 cb\{d) 0 a 2 (<9) 

a 2 (d) 0 a\{d) 0 

0 ci 2 (d) 0 ai(d) 


1 ) 


(3) 


with ai(<9) = — CqV 2 - d 2 , a 2 (<9, cr) = 70 c$d x dy + ^ vd 2 , 

v = cl/c\, and 5 = = 2y/2a^l - T(/// de- 

termines the strength of the coupling between quantum 
fluctuations (set by 1/5) and frustration (set by Ji/J 2 ). 

In the limit A ^ 00 we look for solutions of the type 
cr(x, r) = ctq and iX(x,r) = rag, where do and are 


i 5S' e ff[A,cr 

given by -^ ; 


= 0 . 


= 0 . 


The saddle point equations in the large N limit and for 
the magnetically order phase, where do 7 ^ 0 , then become 


°0 = / ( m 0, 00 ) 

(^) 2 = 


(4) 


We are interested in the quantum phase transition in 
which case I If ? anc ^ SB] 

[(?i(m 0 , (T 0 ) + GZ (m 0 , a 0 ) 


/(m 0 ,cr 0 ) = 1 - g J 
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FIG. 2 : Solutions to Eqs. for the OP ao and the spin 
gap mo/co (inset), at T = 0, as a function of g , for different 
small ratios of J 1 /J 2 . For J\ — 0 (black squares), the OP 
and the spin gap vanish continuously at g c , and the PT is 
of the 2nd order (Heisenberg model). For J\ 7 ^ 0 (red, blue 
and green symbols), however, the OP and the spin gap jump 
discontinuously to zero at g c , indicating a 1st order PT. 



Log(1-g/g c ) 


FIG. 3: Plot of Log(a) x Log( 1 — g/g c ) for different values 
of J 1 /J 2 . A 2nd order PT, at the mean field (large N) level, 
produces a straight line (black squares) with slope [3 — 1 / 2 , 
where a —>• 0 as g —>• g c , as expected for J\ = 0 (Heisenberg 
model). With frustration, Ji ^ 0 (red, blue and green sym¬ 
bols), however, a 7 ^ 0 as g g c (with increasing saturation 
value for a with increasing J 1 /J 2 ) indicating a 1 st order PT. 


h(m 0 ,a 0 ) = y 


/ 


d 3 k 

W) 


3 z 


G- (m 0 , cr 0 ) -Gi(m 0 ,cr 0 ) 


where k z = uj/cq. The Green’s functions are 
Gi(mo, (To) = - —/ —vt, where we have defined 

D^(c j°) = kl + ky ± 70 k x k y + k 2 z (l ± ^). We should 
emphasise now that the unusual coupling between the 
order parameter, ao, and the frequencies, k z = cj/co, in 
(ao) will be responsible for the first order character 
of the quantum phase transition. 

Eqs. @ determine the phase diagram of the model. 
By solving the above set of equations self consistently we 
obtain the behaviour depicted in Fig. [ 2 ] We observe that 
while for the Heisenberg model the OP goes smoothly 
to zero at g c (indicating a 2nd order PT), frustration 
brings the system closer to the QCP and the OP jumps 
discontinuously to zero at g c , indicating a 1st order PT. 
The same is true for the spin gap (see the inset) when 
the transition is approached from the nonmagnetic side. 

To further establish the 1st order nature of the PT, 
we exhibit, in Fig. [3j the dependence of Log(ao) as a 
function of Log(l — g/g c )- For J\ = 0 the behaviour is of 
a straight line with slope given by /3 = 1 / 2 , as expected 
for a mean field behavior (large N) in a 2nd order QPT 
(with the order parameter vanishing continuously at the 
quantum critical point). For J\ 7 ^ 0, however, we observe 
that, although away from the critical point the deviation 
from mean field behaviour is very small, closer to g c the 
deviation is significant and characteristic of a first order 
quantum phase transition, with ao saturating as g g c . 

Let us now provide definitive analytical evidence that 


the transition is indeed 1 st order and not a sharp 2 nd 
order PT. We note that the parameter values obtained 
from the self-consistent equations are such that we can 
write mo /cq as a function of ao fl 6 ] 



(5) 


where 61 (ao) 


(27r) 3 


-D-{* 0 )/ a 2 
e k _ 

DZ(oo) 


e -r>+(<T 0 )/A 2 \ 

6 r 


and 62 (a 0 ) 


-I 



e-°i (J ° )/A2 

(£> 7 ( o - 0 )) 2 


-D+(a 0 )/ A 2 \ 

•cw ■ 


and within such approximation we can rewrite the system 
of self-consistent equations © in terms of a single self- 
consistent variable, namely 


(To = \J /((To5 TOq((To)). 


(6) 


Fig. El shows the plots of Eq. © for Ji/J 2 = 0.2 
and for different values of the coupling constant g. For 
9 = 9i < 9c (green triangles) yi(cr 0 ) = FTFoThiFo)) 
crosses the straight line ^ 2 (^ 0 ) = cro at only one point, 
giving the value of the staggered magnetisation for this 
value of the coupling constant. For gi < g = g 2 < 9c 
(blue diamonds) however we see that 7/1 (ao) and ^ 2 (^ 0 ) 
cross twice. The first (smaller) value of the magnetisa¬ 
tion, however, corresponds to a local maximum of the 
free energy U(a) (inset: unstable fixed point) and shall 
be discarded, while the magnetisation is then determined 
by the second (higher) crossing point solely. By further 
increasing the coupling constant g = go = g c (red circles) 






























4 



a 

FIG. 4: Solution of Eq. © for fixed J1/J2 = 0.2. For g < g c 
(green triangles and blue diamonds) there is only one stable 
equilibrium solution for a (inset: minimum of U(a)). For 
g — g c (red circles), Eq. © produces a nonzero value of cr, 
while for g > g c (pink squares) no solution is found and a 
jumps to zero discontinuously above g c , as in a 1st order PT. 



FIG. 5: Complete phase diagram of the Ji — J 2 model gen¬ 
erated by the NLSM for both Neel and collinear orders, 
0 < a < 1. The Neel order parameter a (filled black squares) 
vanishes continuously indicating a 2nd order PT while the 
collinear order parameter a (filled red circles) jumps to zero 
at the critical point, indicating a 1st order PT. The spin gaps 
mo/u (empty black squares and empty red circles, scaled up 
by a factor 5 for clarity and u = \^2SaJi) cross linearly at the 
classical border at a m 0.5, in agreement with DMRG m 

we find a single critical solution to Eq. © giving a finite, 
nonzero and sizeable value for the staggered magnetiza¬ 
tion, which, however, ceasses to exist for g = g 4 > g c 
(pink squares). The fact that the sublattice magnetisa¬ 
tion jumps to zero discontinuously for g > g c indicates 
the 1st order nature of the QPT. 

It is important to emphasise that the magnon disper¬ 
sion along the collinear directions k x = k y = k is given 


by wjvLSM(fc) = c_a/2 - To |fc|, with c _ 2 = c 0 2 - c, 2 /2 
m he and thus acquires an imaginary part beyond 
the border at a < 0.5 when 70 > 2, showing that the 
magnetic excitations of the collinear state move from 
q m (77 0) and/or q = ( 0 , 7 r) towards the one of the Neel 
state at q = (77 7 r), as expected. The complete phase 
diagram obtained within the NLSM formalism, for the 
whole range 0 < a < 1 is given in Fig. [5] 

We have obtained the complete phase diagram of the 
antiferromagnetic J 1 -J 2 Heisenberg model within the 
framework of the O(N) nonlinear sigma model. We have 
found that the two magnetically ordered phases, Neel or¬ 
der for a < 0.4, and collinear order for a > 0.6, are 
separated by a nonmagnetic region at 0.4 < a < 0.6 
where a gapped spin liquid is found. The transition at 
a = 0.4 is of the second order while the one at a = 0.6 
is of the first order and the spin gaps cross linearly at 
a = 0.5. Our results are exact at N 00 and agree with 
numerical results from different methods. 
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